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We present observations of the insulator-superconductor 
transition in ultrathin films of Bi on amorphous quartz, quartz 
coated with Ge, and for the first time, solid xenon condensed 
on quartz. The relative permeability tr ranges from f .5 for Xe 
to 15 for Ge. Though we find screening effects as expected, 
the TS transition is robust, and unmodified by the substrate. 
The resistance separatrix is found to be close to h/4e^ and 
the crossover thickness close to 25 A for all substrates. I — V 
studies and Aslamazov-Larkin analyses indicate superconduc- 
tivity is inhomogeneous. The transition can be understood in 
terms of a percolation model. 

PACS Numbers : 73.50.-h, 74.40. +k, 74.80.Bj 

The insulator-superconductor (TS) transition has been 
extensively investigated over t|he.last decade, in a variety 
of syiStems such as thin filmailH single Josephson junc- 
tional, arrayso, one dimensional wireaJ. Values of limit- 
ing resistance close to the quantum resistance for pairs 
in one dimensional wiresO, and two-dimensional, nomi- 
nally homogeneous ultra-thin filmsEJ are reported. Dif- 
fering values of thn^limiting resistance at the transition 
have been observedo in different systems, and attributed 
to structure, i.e. homogeneous and granular films are 
expected to behave differently. 

A-l phase-only picture, first proposed|-|by Ramakrish- 
nanQ, and further elaborated by FisherEl has been con- 
sidered appropriate for such systems. A scaling theory 
of the I-S transition has been developed.H This theory 
predicts that the critical resistance will be universal, i.e. 
independent of all microscopic details, if the system is in- 
variant under the interchange of the roles of chajEge and 
flux. Thus the TS transition may be self-dualEl. The 
precise value of the critical resistance appears to depend 
on the nature of the interaction between charges, being 
equal to h/Ae^ only when the interaction is logarithmic 
in the separation, one of the conditions which must exist 
for self-duality. The interaction between charges can be 
logarithmic in their separation if th^ 2D films have a suf- 
ficiently high dielectric constant JlJ This condition may 
be met for semiconductor and semimetal films. 

This has to be reconciled with experimental observa- 
tions, which show that superconductivity in ultrathin 
films canJac enhanced by the proximity of a metal or 
diclcctricEJ This is attributed to partial screening of the 
Coulomb interaction between conduction electrons in the 
films. Screening also affects the properties, of insulat- 
ing films, by changing the localization length.ll3llj Glover 
studied the effect of the substrate dielectric constant on 
transition temperature of thin films, over thirty years 
ago, inconclusively. Ge or Sb has been |extfiiisively used 
as underlayers in ultrathin film studieil'E£lE3. A natu- 



ral question arises at this point about the effect of the 
substrate/underlayer on the TS transition. 

For a film on a under layer, the general expression for 
the interaction potential between charges in the film is 
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where d is the screening distance (the distance between 
the midplane of the film and the midplane of the under- 
layer) and r the separation between the charges in the 
film. If the separation of two charges is small compared 
with the distance between a charge and its image (2(i), 
the underlayer makes little difference, and the interac- 
tion remains monopolar. At large distances, however the 
charge and its image behave as a dipole and interaction 
falls off more rapidly with distance. For the films studied 
here, 2d would be close to the film thickness. Although 
screening by charge carriers in the film is considered in 
the usual BCS treatment, the rclatiiie permeability for 
a highly disordered film is unknown.EII It is expected to 
lie between the metallic and insulating limits of — oo 
and €r — 10. This issue becomes important in the TS 
transition region. U{q), the term in the BCS interaction 
Hamiltonian becomes 



47rp^ 

q + kj 



(2) 



where kg is the inverse screening length. In principle a 
film with a transition temperature of K defines Rc- As 
we show, screening influences the transition temperature 
of superconducting films, and the normal state resistance 
of all films. Therefore it should affect the TS transition. 

It has also been proposed that the TS transition in the 
limit T = is a combined effect of pairing and local- 
ization in 2D systems.lla Experiments have refuted this 
conjecture. citekagawa How disorder, cluster size, and the 
Coulomb interaction influence the behavior of Cooper 
pairs is still not very clear. It has been conjectured that 
the Coulomb interaction majj-^ot influence the critical 
normal state sheet resistance.tl However, this is refuted 
by observations of the M-I transition in siliccp, MOS- 
FETS and other systems in zero magnetic field.c3 It has 
been shown that interactions may increase or decrease 
the conductance of a disordered 2D electron system - 
weak interactions increase the dc conductance in the lo- 
calized regime while-they decrease the conductance in 
the diffusive regime. Ej Strong interactions were always 
found to decrease the conductance. These considerations 
motivated our experiments. 

In superconductors. Coulomb interactions suppress the 
fluctuations in the number of electrons, and increase the 
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fluctuations of the phase of the superconducting order 
parameter. This affects the Josephson coupling energy 
Ej in inhomogeneous systems, where the competition 
between the charging energy Ec and Ej drives the tran- 
sition, citekagawa Typical systems exhibiting such behav- 
ior are granular films, and junction arrays. Ej may 
be estimated using the Ambegaokar-Baratoff equation, 
Ej = 7r^A/4e^i?Ar, where A is the BCS gap for Bi, Rn 
is the normal state resistance of the film. Ec can be esti- 
mated as Ec = q'^/ineoKd where k = er{d/2s + l), where 
er is the relative permeability of the substrate, d is the 
size of the grains/clusters, and s is the spacing. Since Rn, 
and K depend on the substrate/underlayer, one expects 
the I-S transition to occur at different Rn on different 
substrates. Since the effects of disorder, cluster size, and 
interactions on Cooper pairs are unclear, exlpicit predic- 
tions of Rn for the I-S transition on various substrates 
cannot be made. 

In this work, we report studies of the I-S transition in 
quench condensed Bismuth films, as a function of disor- 
der (or film thickness) on a variety of substrates - amor- 
phous quartz, quartz coated with Ge, and solid xenon 
condensed on quartz. The relative permeability e,- ranges 
from 1.5 for Xe to 15 for Ge. Despite screening effects 
the I-S transition is robust. Studies on other substrates 
such as single crystal sapphire, and sapphire coated with 
Ge, Xe yield similar results. The experiments were done 
in a custom UHV cryostat equipped with reflection elec- 
tron diffraction (RHEED) , capable of a hydrocarbon free 
vacuum of ~ 5x10^^" Torr, described elsewherecj. The 
Ge underlayers were deposited on the substrates in a sep- 
arate UHV system at room temperature. The Xe under- 
layers were grown in-situ. Our experimental setup re- 
sembles that of Ref.[l-2], although we cannot attain such 
low temperatures. We quench condense our ultra-thin 
films in the temperature range 1.8 K to 15 K. The sub- 
strate temperature influences the disorder in the film, 
and thereby its properties. Here we report results on 
films that were quench condensed at 15-iK, which facili- 
tates comparison with published resultsa. 

RHEED studies show that the Bi is almost amorphous. 
It is difficult to unambiguously distinguish between amor- 
phous and nanocrystalline at such low temperatures, on 
poorly conducting films due to charging effects. Based on 
the Scherer formula for the peak broadening, we estimate 
that films thicker than 10 A are ccpiposed of clusters that 
vary in size from 25 A to 100 AE3. Since the informa- 
tion obtained is in reciprocal space, it is difficult to com- 
ment on the real space surface morphology. Our RHEED 
observations are consistent with previous results. E3 Su- 
perconductivity in granular systems of rhombohedral Bi 
clusters has also been reported.tj 

Figure 1 shows the evolution of the temperature depen- 
dence of the sheet resistance R(T) with thickness for (a) 
Bi films on Ge (10 A thick), which has been deposited on 
amorphous quartz, (b) Bi films on quartz, and (c) Bi films 
on solid xenon condensed on amorphous quartz. A tran- 
sition from insulating type behavior, to superconducting 



behavior as the thickness of the films is increased is clear. 
This type of zero field transition is considered a zero tem- 
perature quantum phase transition, controlled either by 
disorder, carrier concentration, or thickness. The normal 
state resistance at an arbitrarily high temperature Rn 
has traditionally been used to parameterize the transi- 
tion, although it may be weakly temperature dependent 
above the superconducting transition temperature, and 
becomes ill defined as the I-S transition is approached. 
The value of the normal state resistance of a film on the 
boundary between superconducting and insulating be- 
havior has been referred to aS||tlie resistance separatrix, 
and has been denoted by Rjdffl. We obtain Rq as an 
algebraic average of the sheet resistances of the last insu- 
lating and the first superconducting films, measured at a 
relatively high temperature (10 K). Rq is close to h/Ae^ 
for all three sets of data. This observation indicates that 
the value of Ro is substrate independent, and possibly 
experiment independent. 

The transition temperature for a film of a given thick- 
ness is higher on substrates of higher relative permeabil- 
ity, and the normal state resistance is also higher. A 65 
A film on Xe, quartz and Ge has T^s of 3.8, 4.2 and 
4.42 K respectively. These results are consistent with 
published dataJiJ This shows that the conductance of a 
disordered film-depends on interplay between interaction 
and disorderO We caution that since our experiments 
are limited to 1.8 K, we cannot rule out the possibility 
that a film which appears to be the last film on the in- 
sulating side, maji, at lower temperature turn out to be 
superconducting.ll3. A plot of conductance vs thickness 
at different temperatures gives the crossover thickness, 
which is shown in Fig. 2 for films on Xe. We find the 
crossover thickness to be between 25 to 28A for all sub- 
strates. 

The insulating films follow a behavior that is consistent 
with the results reportcd.Ed We find that the conductiv- 
ities of all our insulating Bi films, on all substrates to 
be: 
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where x changes with film thickness. For the thianest 
films, X is close to 0.75, due to collective hoppin^3. As 
the thickness is hivthex. increased, x reduces to 0.5, the 
Efros-Shklovskii formE3, which describes hopping modi- 
fied by Coulomb interactions, and finally for the thjckest 
insulating films x is close to 0.33, the Mott valueEZl. 

The form of the R(T) for these films may lead us to 
the conclusion that these films are nominally homoge- 
neous. Such a conclusion is incorrect as we shassj below. 
Although STM studies of surface morphologycJ report 
that even at 75% coverage the films are not conduct- 
ing, hence no evidence for percolative behavior, we find 
that supepcpnductivity in our filmsJ^i indeed percolative 
in naturcEd Aslamazov and Larkina considered the pos- 
sibility of fluctuations causing superconductivity. The 
total conductivity is given by ct = cr at + cr2_D' where un 
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is the normal state dc conductivity, and a the paracon- 
ductivity. Its temperature dependance is similar to that 
of the magnetic susceptibility at T>Tc. They derived 
the result. 
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mean field transition temperature. t/Rq is a constant 
(<7AL=e^/16?i) for all materials. We have evaluated t/R,^ 
idexp) for various films. A systematic dependence of gexp 
on the thickness d is shown in Fig. 3. This parameter 
deviates from gAL for thinner films. It approaches the AL 
value (gAL) as the thickness is increased. It is assumed 
that theory predicts gexp = gAL for all films, independent 
of microstructure. Both the normal state conductance 
and paraconductance depend on sample shape. GlovertJ 
has shown that as the microstructure deviates from a 
"uniform rectangular slab", gexp exceeds the AL value. 
The thinner the film, higher the disorder, larger are the 
deviations from a slab geometry, and larger the deviation 
of gexp from the AL value. Hence, films close to the 
transition are inhomogeneous. 

Further evidence is in the form of hysteretic I-V curves, 
shown in Fig. 4. These curves can be understood in 
terms of a resistively and capacitively shunted random 
Josephson junction array model. These I-V^ have been 
discussed in detail in a separate publicationEHl. The phys- 
ical picture is of superconducting islands connected by 
thin regions of normal metal, which act as the weak links 
between the islands, thus forming a random JJ array. 
The distribution of grain/cluster, each with a finite num- 
ber of electrons, requires-a finite N BCS treatment. This 
results in a spread in Tj2l, so that at some temperature 
all regions of the film arc not superconducting. Normal 
regions exist, and act as weak links. Areal maps of quasi- 
particle conductance or superconducting order parameter 
amplitude can resolve such regions. j— . 

Feigel'man and LarkirH and Spivak et a/O have dis- 
cussed the quantum superconductor-metal transition in 
a 2D proximity coupled array. Thickness larger than the 
coherence length was considered by the former, with the 
opposite extreme analyzed by the latter. They predict 
intervention of a normal phase. Limitations on the ac- 
cessible temperatures prevent our study confirming this. 
From the Hc2(T) curve, we have determined the coher- 
ence lengths (^c) to be close to the film thickness. This 
regime merits further theoretical study. We have eval- 
uated the localization lengths (^/), and we find that at 
the I-S transition, ^; >> ^c, consistent with the results 
ofRef. [17]. 

The effect of strong disorder on the superconducting 
order parameter amplitude has been investigated within 
the Bogoliubov-de Gennes frameworkl3. The local pair- 
ing amplitude develops a broad distribution with signifi- 
cant weight near zero, as disorder increases. The density 



of states showed a finite spectral gap. Persistence of the 
gap was found to arise due to the breakup of the system 
into superconducting islands. Superfiuid density and off- 
diagonal correlations showed a substantial reduction at 
high disorder. Incorporation of phase fiuctuations lead 
to a non-superconducting state. Our data are consistent 
with such a picture. 

Our results suggest that a percolation type model, pro- 
posed by MeirtJ to explain the M-I transition in 2D sys- 
tems is most appropriate and merits investigation. Meir 
suggested the following physical picture for the M-I tran- 
sition in 2D - potential fluctuations due to disorder de- 
fine density puddles of size or larger, within which the 
electron wave function totally dephases. Locally between 
these puddles, transport is via quantum tunneling. Sup- 
port for this idea comes from fits of the conductance at 
a high temperature (10 K) for the insulating films which 
follows a power law, with an exponent of 1.33 charac- 
teristic of 2D percolation, as shown in Fig. 5. As an 
interesting observation, we also show that the supercon- 
ducting films also fit this same power law, albeit with 
a different critical thickness. The critical thickness for 
the insulating films is the thickness for onset of electri- 
cal conduction. We have no explanation as to why the 
critical thickness for normal state conduction depends 
on whether or not the film superconducts at lower tem- 
perature, but present this as an interesting observation. 
Keeping in mind the model discussed above, it is logi- 
cal to associate the puddles with superconducting regions 
and the quantum point contacts with the normal regions. 
The competition between the Ec and Ej probably drives 
the transition. However, it is intriguing that the transi- 
tion is robust at /i/4e^ on all substrates, despite Ec and 
Ej being depepdcpt on e^. Percolation issues have been 
studied earlier.E3"c3 

The conditions for 2D coulomb interactions may-have 
been met in this and other experiments.EJ Keldyshli3 has 
shown that this requires the dielectric constant of the 
film to be much greater than that of the substrate. In 
this study, substrate dielectric constants range from 1.5 
for solid Xe to 15 for Ge underlayers. Hall measure- 
ments were done for the first superconducting film near 
the transition, and this yielded a high carrier concentra- 
tion of the order of 10^^cm~L^ which is consistent with 
that reported in the literatureEZl, suggesting a high dielec- 
tric constant for the films. We have not done Hall mea- 
surements for films on the insulating side, since interpre- 
tation of data in the hopping regime is complicated. Vor- 
tices exist in these fUnis, since they are random Joseph- 
son junction arrays,cll due to the discrete nature of the 
array, which requires that the current take convoluted 
paths, consistent with a transition governed by percola- 
tion. Duality in JJ arrays has been observed before.Q 

In conclusion, we have observed robust I-S transition 
in ultrathin Bi on several substrates. We have presented, 
for the first time, the I-S transition on solid inert gas 
underlayers. Although the substrate infiuences the Tc as 
well as Rn of the films, the transition itself is indepen- 
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dent of the interaction potential between charges. An 
AL analysis indicates inhomogenous films, contrary to 
the R(T). I-V's indicate a percolation type transition. A 
model similar to that proposed by Meir for the M-I tran- 
sition in 2D merits further investigation. Further studies 
at lower temperatures are needed. 
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FIGURE CAPTIONS 
Figure 1 Evolution of Rq vs T of Bi on (a) 10 A Ge 
on amorphous quartz (b) bare quartz and (c) solid Xe 
condensed on amorphous quartz. The substrate temper- 
ature during deposition in all cases was maintained at 15 
K. 

Figure 2 Plots of resistance vs. thickness at different 
temperatures showing the crossover thickness of about 
26A for Bi films on Gc. The crossover occurs within 
ibSA of this value for all substrates. 
Figure 3 The parameter g^^p which characterizes the 
amplitude of the fluctuations as a function of film thick- 
ness. It asymptotically approaches the Aslamazov- 
Larkin value. 

Figure 4 I-V of the superconducting films on quartz. 
The hysteresis loops show the dissipation due to resis- 
tances shunting the junctions of the random Josephson 
junction array. The first film is 40A, and the thickness 
increment is lOA. 

Figure 5 Fits of the conductance of the insulating films 
(taken at 10 K) and the normal state conductance of 
the superconducting films (also taken at 10 K) to 2D 
percolation function, for Bi films on Ge. Only the critical 
thicknes in each regime is different, as discussed in the 
text. 
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